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Abstract 

We analyze the influence of the density dependence of the symmetry energy on the average ex- 
citation energy of the isoscalar giant monopole resonance (GMR) in stable and exotic neutron-rich 
nuclei by applying the relativistic extended Thomas- Fermi method in scaling and constrained cal- 
culations. For the effective nuclear interaction, we employ the relativistic mean field model supple- 
mented by an isoscalar-isovector meson coupling that allows one to modify the density dependence 
of the symmetry energy without compromising the success of the model for binding energies and 
charge radii. The semiclassical estimates of the average energy of the GMR are known to be in 
good agreement with the results obtained in full RPA calculations. The present analysis is per- 
formed along the Pb and Zr isotopic chains. In the scaling calculations, the excitation energy is 
larger when the symmetry energy is softer. The same happens in the constrained calculations for 
nuclei with small and moderate neutron excess. However, for nuclei of large isospin the constrained 
excitation energy becomes smaller in models having a soft symmetry energy. This effect is mainly 
due to the presence of loosely-bound outer neutrons in these isotopes. A sharp increase of the 
estimated width of the resonance is found in largely neutron-rich isotopes, even for heavy nuclei, 
which is enhanced when the symmetry energy of the model is soft. The results indicate that at 
large neutron numbers the structure of the low-energy region of the GMR strength distribution 
changes considerably with the density dependence of the nuclear symmetry energy, which may be 
worthy of further characterization in RPA calculations of the response function. 
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I. INTRODUCTION 



The departure of the equation of state (EOS) of neutron-rich matter (A^ > Z) from 
the symmetric hmit (A^ = Z) is governed by the nuclear symmetry energy. This makes the 
knowledge of the symmetry energy important for a variety of areas of nuclear physics (masses, 
neutron densities, heavy ion collisions, etc.) and of astrophysics (supernovae, neutron stars, 
nucleosynthesis, etc.) jll-H]. However, the available experimental data on nuclei, which 
correspond mostly to isotopes with small and moderate neutron-proton asymmetries, do not 
constrain very precisely the value of the symmetry energy at the saturation density. Rather, 
it seems to be better constrained at a lower density around 0.10 fm~^ [Z lol|, corresponding 
to an average between the bulk and surface symmetry contributions in the semi-empirical 



mass formula 11-13 



^ ^ All in all, the density dependence of the symmetry energy is a 

basic nuclear property that still is not sufficiently well known [l|, 0|. This situation has 
prompted a recent ongoing effort to characterize more closely the symmetry energy in the 
regime of subsaturation densities, relevant for finite nuclei, from the phenomenology of 
different nuclear observables [H, 0, ll|"19|. A proper understanding of the symmetry energy 
is particularly desirable in order to predict ground-state and excitation properties of atomic 
nuclei far from the naturally-occurring region of stability (20-24 . 



The study of collective excitations of nuclei near and away from the stability line is 



For instance, the isoscalar 



a salient source of information for nuclear structure (25-291 
giant monopole resonance (GMR) constrains the compression modulus Kq of the EOS of 
symmetric nuclear matter. This resonance has been accurately measured in heavy and 
medium mass nuclei through recently improved a-particle scattering experiments 27|, |28 



By comparing the Ip-lh RPA results predicted by Gogny [3^], Skyrme |3lH33| and relativistic 
mean field (RMF) interactions (34, 35| with the excitation energy measured in ^osp]^ (14.17± 
0.28 MeV [27|), it is concluded that the value of the nuclear incompressibility Kq must lie 
in the range of 200-280 MeV. Performing measurements of the GMR in exotic nuclei may 
deliver other valuable informations and enhance our knowledge of the physics of isospin-rich 
and loosely-bound nuclear systems 2^, 2^, 32, 36, 37]. However, it represents an important 
experimental challenge due to e.g. the low intensities in the radioactive ion beams. With 
the advent of new technologies it is hoped that measurements of the breathing mode in 



short-lived unstable nuclei will become feasible in the future [26|. Note that promising first 
successes have already been demonstrated [29j]. Although the value of the excitation energy 
of the GMR is driven to a large extent by the compression modulus of symmetric matter, 
it is known to be infiuenced also by the density dependence of the symmetry energy in a 
neutron-rich nucleus such as ^osp^^ i 



38|-|42|. This infiuence is expected to be amplified 



in isotopes with large neutron numbers, but how and to which extent has not actually been 
assessed. 

Our purpose in this paper is to investigate the effects of the density dependence of the 
symmetry energy on the average excitation energy of the GMR in isospin-rich nuclei, of 
heavy and medium mass, when moving away from the stability valley. We ought to mention, 
however, that our aim here is not pursuing to constrain the nuclear EOS from the present 
calculations, but rather to analyze the effects that arise from the existing uncertainty in 
the knowledge of the density dependence of the symmetry energy. We advance that the 
calculations will show that the impact on the excitation energy of the GMR is somehow 
moderate for isotopes not far from stability, while the more notable effects are found for 
nuclei with large neutron numbers, thereby (unfortunately) rendering them more difficult to 
see in experiment. To describe the effective nuclear interaction, we found convenient to use 
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the model of jS-lO] that consists of a standard RMF Lagrangian plus an isoscalar-isovector 
nonlinear interaction between the u and p meson fields. This coupling provides an efficient 
way to change in a controlled manner the density dependence of the symmetry energy 
without modifying the EOS of symmetric matter. In particular, the compression modulus 
Kq is not modified. When the model is applied to finite nuclei, the variation of the strength 
of the isoscalar-isovector coupling changes the properties of the neutron density distribution 
(for example, the neutron rms radius), which have large experimental uncertainties, but 
the predicted ground-state binding energies and charge rms radii, which are accurately 
known in experiment, remain basically unaltered. Thus, this RMF model allows one to 
modify the density content of the symmetry energy while preserving the agreement with the 
accurate grou nd-state information on which the model has been calibrated previously (see 



mm 41 



The consistent microscopic framework for the calculation of the properties of the GMR is 
furnished by the Hartree-Fock plus RPA theory, which corresponds to the small-amplitude 
limit of the time-dependent Hartree-Fock method. The generalization of this theory to 
deal with pairing in open-shell nuclei is known as quasiparticle RPA (QRPA). The QRPA 
provides the full information on the response function and centroid energies of the GMR. 
These microscopic calculations are manageable but rather demanding if they are fully self- 
consistent, which is accentuated when they are performed using the canonical basis or in 
isotopes with large neutron numbers as it has been recently discussed in [s^l . We will carry 
out our analysis by studying the estimates of the average excitation energy of the GMR that 
can be obtained without performing explicit RPA calculations but which carry meaningful 
physical information on the RPA strength. The typical example of this procedure is the 
sum-rule approach to the RPA developed in the nonrelativistic framework using Skyrme 



forces 144 - 46 



The sum-rule techniques provide average properties of the GMR such as the 
centroid energy and the resonance width and have been successfully applied to study the 

In the sum-rule approach 



33,^ 



GMR in different nuclei and isotopic chains [3C 
one obtains the cubic and inverse energy-weighted moments of the RPA strength function 
by means of, respectively, scaling and constrained calculations of the uncorrelated Hartree- 
Fock nuclear ground state. Note that these ground-state scaling and constrained values 
are the exact RPA moments j44-46|. Because highly-accurate ground-state calculations are 



simpler, the sum-rule techniques also have been used in the literature as a means to check 



the self-consistency of explicit RPA calculations of the strength function [32|, |37|, |47|, |48 



The validity of the sum-rule theorems when pairing is active (QRPA) has been proven 



in recent publications [37|, |49| . Capelli et al. [37[ have made a careful study of constrained 
Hartree-Fock (CHF), constrained Hartree-Fock-Bogoliubov (CHFB) and QRPA calculations 
of the excitation energy of the GMR in the light isotopic chains of Ca and Ni. The calcula- 
tions have reached up to the theoretical neutron drip line of the Skyrme SkM* force, i.e., up 
to the isotopes ^^Ca and ^^Ni. Admitting that these isotopic chains are known to terminate 
earlier in experiment than in the predictions of SkM* and of most of the nuclear effective 
forces, and knowing also that it will be difficult to measure the GMR in even less exotic 
isotopes, the calculations spanning the whole theoretical isotopic chains are very useful to 
test the predictions of nuclear forces under extreme isospin conditions (which can take place 
in astrophysical scenarios [50|), as well as to test the different many-body frameworks em- 
ployed to deal with the problem. It may be noted also that Khan 5l| has recently concluded 
that knowing the GMR in a single isotope such as ^osp]^ ^[y^y ^^^^ suffice to extract Kq of 
nuclear matter and that the whole isotopic chain should be measured in order to have a 
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grasp on the various effects acting on tlie GMR [51 



From tlie results of [37(| for Ca and Ni, it is seen that the inclusion of pairing (CHFB 
and QRPA calculations) does not modify significantly the behaviour shown by the GMR 
excitation energies obtained in the CHF calculations along the isotopic chains. In the case 
of the Ni chain, the atomic number not being as small as in Ca, the CHFB and CHF values 
agree fairly well for most isotopes from ^^Ni up to the neutron drip line included, as noticed 



in figure 5 of [37[. Of course, superfluidity is to be taken into account if one wants to 
describe in detail the excitation spectra or if precise values of the excitation energy within a 
few hundred keV in specific nuclei are needed for comparison with accurate measurements 



371 . |51| . |52| . It has been pointed out in [37[ that the QRPA calculations for isotopes having 
large neutron numbers become much time consuming and that ensuring proper numerics 
in such isospin-rich nuclei becomes somewhat difficult due to the weakly bound neutrons. 
Indeed, the authors of state that in the case of large neutron numbers very accurate 
QRPA calculations may become almost prohibitive in keeping with the size of the basis. 

The GMR is a collective excitation whose average properties in general vary smoothly 
with the mass number A and are rather insensitive to shell effects, at least for not too light 
isotopes. Thus, a description of the average properties of this resonance through semiclassi- 
cal approaches like the extended Thomas-Fermi (ETF) method can be, for some purposes, 
a practical and interesting alternative. It has been applied in different studies in the non- 
relativistic Skyrme-Hartree-Fock framework [i^, 47, The results have shown that 



the average excitation energies of the GMR obtained from scaling and constrained quantal 
HF calculations are nicely described by the corresponding semiclassical ETF counterparts. 
The ETF predictions of GMR excitation energies become increasingly closer to the quantal 
predictions with increasing atomic number, as it can be expected in semiclassical methods. 
Indeed, in the calculations along the Pb isotopic chain with the Skyrme force shown in 
figure 12 of [13], the ETF scaling and constrained excitation energies differ from the cor- 
responding quantal values by less than 0.4 MeV and they reproduce quite accurately the 
change of the quantal results in going from the proton to the neutron drip line. Thus, from 
the experience in nonrelativistic calculations, we believe that the semiclassical scaling and 
constrained estimates of the average excitation energy will be a useful and realistic procedure 
to have a quick representation of the general trends of the GMR along isotopic chains also 
in the relativistic mean field framework. And in heavy systems such as the Pb isotopes, we 
expect that the semiclassical calculations of the GMR average excitation energy will yield 
reliable predictions practically on the quantitative level. Certainly, for a precise microscopic 
description, the more demanding QRPA calculations would have to be performed. 

The semiclassical relativistic extended Thomas- Fermi (RETF) approach [56-60] has re- 
cently been applied in scaling and constrained calculations of the GMR energy in stable 
isotopes in |6l|, |62| . The results turn out to be in good agreement with the values extracted 



from full relativistic RPA calculations l6lL 621. Let us mention as a side comment that the 



RETF scaling procedure has been helpful in obtaining the virial theorem for the relativistic 
mean field model j6l| and that it has allowed self-consistent calculations of the surface con- 
tribution -ft'surf to the nuclear incompressibility in the relativistic theory j63|. In this paper 
we will use the relativistic nuclear mean field model with an isoscalar-isovector coupling 
mentioned previously together with the RETF approach for describing the average excita- 
tion energy of the GMR in the Pb and Zr isotopic chains, through scaling and constrained 
calculations. The present results are expected to provide a realistic indication on how the 
density dependence of the symmetry energy impacts on the behaviour of the excitation en- 
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ergies and widths of the GMR along isotopic chains of heavy and medium mass. In the next 
section the basic theory is presented, in the third section the results are discussed and the 
summary and conclusions are laid in the fourth section. Some details about the semiclassical 
RETF calculations are given in the appendices. 



II. FORMALISM 



A detailed description of the scaling and constrained calculations of the excitation energy 
of the GMR in the relativistic framework using the RETF approach has been presented 
elsewhere 61, Here, we shall restrict ourselves to outline the more relevant expressions, 
incorporating the new contributions that arise from the mixed u-p meson interaction. 



A. ETF energy density of the relativistic mean field model 

The basic model for the present study is the relativistic mean field theory of nucleons and 
mesons. As in the nonrelativistic problem, the theory can be formulated in the Thomas- 
Fermi approach. When corrections of order /i^ (in gradients of the nucleon densities and of 
the Dirac effective mass) to the pure Thomas-Fermi contribution are included in the energy 



density, it is known as the relativistic extended Thomas-Fermi approach |56l-l60 



The local energy density that we employ in the present work can be written as follows: 

- ^ [i^Wf + mlW^] - ^ [{VBf + mlB^] - AyB^W^ - i (V^)^ , (1) 

where p = pp + pn is the baryon density, ps = \{pp — Pn) is the isovector density, and the term 
£ will be described immediately below. The energy density ([1]) contains scalar $ = gs(f)o{f), 
vector W = gvVo^r) and vector-isovector B = gpbQ^r) meson fields, as well as the Coulomb 
field A = 6^40(7"). These fields represent, respectively, a- meson, w-meson, p-meson and 
photon exchange. The cubic and quartic A<l>^ self-interactions of the scalar field soften 
the EOS of symmetric nuclear matter around the saturation point. In addition, the model is 
supplemented by a nonlinear mixed coupling between the w-meson and the p-meson fields. 
This AyB'^W^ interaction was introduced in to soften the density dependence of the 
symmetry energy and of the EOS of neutron matter, which are two quantities that are 
predicted to be stiff by the majority of relativistic mean field models. 

In the quantum approach, the contribution 8 to the RMF energy density in equation ([1]) 
is given by £^ = ^ + 7o("^ — ^)]fu- Here, the ipu are single-particle Dirac 

wave functions, a and 70 are the usual Dirac matrices and m denotes the rest mass of the 
nucleons. The RETF representation of S consists of a pure Thomas-Fermi part Sq plus a 
part S2, which is of order h"^ compared to Sq. The nucleon variables of the RETF energy 
density are the neutron and proton densities (p„ and Pp) instead of the single-particle wave 
functions (f^. Thereby, £^ is a functional of the nucleon densities and of the Dirac effective 
mass m* = m — $ in the RETF approach. The gradient corrections contained in the S2 term 
provide an improved description of the nuclear surface with respect to the Thomas-Fermi 
treatment with only the Sq contribution. The reader can find the full RETF expression 
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of the functional £ in appendix A. We also provide there the variational Euler-Lagrange 
equations derived from equation ([1]) in the RETF approach for the nucleon densities and for 
the meson and photon fields. 

The energy density "H is to be integrated over space to compute the total energy of the 
nucleus. By means of partial integrations, and on account of the variational meson field 
equations given in appendix A, it is possible to write the relativistic energy density of a 
finite nucleus in the following form: 

n = 8 + Upf + ^<l>^ + + 1 + l^^^ + p^^b2^2 ^ 1^^^ _ (2) 
Here, we have introduced the effective scalar density defined as 

Pf = Ps- - (3) 

The form given in equation ([2]) for H turns out to be more practical to perform the scaling of 
the equations in order to compute the excitation energy of the GMR in the scaling formalism. 



B. Energy of the GMR in the scaling approach 

The excitation energy of the isoscalar giant monopole resonance of a finite nucleus of 
mass number A is customarily written as 



^" = Visr' 

where Ka is called the compression modulus or incompressibility of the finite nucleus, and 
Bm is called the mass or inertia parameter of the monopole oscillation. In the scaling 
approach, the incompressibility is obtained from the restoring force of the monopole 
vibration modeled through a scaling transformation of the nucleon densities: 

Pqair) = a^pg{ar). (5) 

Here, a denotes the scaling parameter and g = n, p refers to neutrons or protons. The meson 
fields do not scale as simple powers of a because of the finite-range character of the meson 
interactions; we collect in appendix B the variational field equations for the scaled meson 



fields of the RMF model with the mixed u-p interaction. Following [6l|, |62|, it is helpful to 
write the scaled Dirac effective mass jtj,* (r) = m — $Q,(r) in the form 

■m*^{r) = a m* (ar) (6) 

because the energy density £ and the scalar density ps of the RETF model scale, respec- 
tively, as Sa{r) = a"^ £[pg{ar) , rh* {ar)] = a'^£{ar) and as Psa{f) = Ps[Pq{c(i^) , ^* {c(f)] = 
a^Ps(ar). The tilded quantities S and ps have the same expressions as £ and ps given in 
appendix A if one replaces m* by m* in them. Taking into account the above, the scaled 
form of the energy density of the present model becomes 



l-inir) = a 



aS + ^<!>^pf + ^^^^<^l + ^^^^^^^ 



(7) 
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with the definition p'f = 

depend on the variable ar. 

The calculation of the first and second derivatives of the scaled e nerg y E[a) 



2! a3 



3! a3 



All densities and fields on the r.h.s. of 



with respect to a proceeds along the same lines described in |6ll . |6 
the scaling incompressibility is 



/ drHair) 



The final result for 



K 



seal 



1 

I 



d{ar) 



Ho 



1 
A 



- / dr 



dps , „ 

-m— h 3 

oa 



9s 9i 9p 



2^$ + -<l>^ 

9l 2 



da 



mt 



+ 2^W 
9v 



dWa 

da 



3! 

.dB, 



9p da 



(8) 



a=l 



where {dps/ da)a=i = —{dps/dm*) [m*+{d(^ a I da) o.=i\. It is interesting to note that all of the 
nucleonic contributions to K^^^"^ are summarized in the term with the scalar density and that 
the massless photon field and the quartic meson couplings such as and KyB^W^ do not 
provide explicit contributions. The expression ([H]) requires the calculation of the derivatives 
of the meson fields with respect to the scaling parameter a at equilibrium (a = 1). These 
derivatives can be computed by differentiation of the scaled meson field equations with 
respect to a and by solving them at a = 1, as we describe in appendix B. 

We will use the notation E^^^ to refer to the excitation energy of the GMR calculated in 
the scaling approach, namely, 



E 



seal 
M 



(^j^seal 



B 



M 



The mass parameter of the monopole resonance is given by the expression [6^, IGJ, l65 



B 



M 



dr r^U . 



We note that in the nonrelativistic limit it reads 



dr r'^mp = m A (r^) , 



(9) 



(10) 



(11) 



where m is the nucleon rest mass. Because of the large contribution of the nucleon rest mass 
term mp to the energy density "H in the integrand of equation (ITU]) , the value of the excitation 
energy of the resonance is httle modified by using either Bm 



or 5^ in the calculations. 



C. Energy of the GMR in the constrained approach 



One also can study the avera ge e xcitation energy of the giant monopole resonance through 
constrained calculations 44, 6^, 6^6^. For that purpose, one solves the problem associated 
with the constrained functional 



j dr[H — Tjr'^p] = E{7]) ~ V J drr' 



(12) 
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The densities, fields and energy obtained from tlie solution of the variational equations 
deduced from equation (fT2|) depend on the value of the parameter t] at which one performs the 
calculation. By expanding -^(r/) in a harmonic approximation about its minimum, located 
at the ground-state rms radius Rq (corresponding to = 0), one obtains the compression 
modulus of the finite nucleus in the constrained approach: 

^r'^i«?f^) . (13) 



Here, R"^ = (r^)^ is the square mass radius of the nucleus computed with the density p of 
equation ( !T2|) . One can show that equation (fT3|) can be rewritten in the following equivalent 
forms: 

Kr-^-^KC-^Y -4ARi('-W)\ (14) 



We have used f|T4|) to check the numerical accuracy of our calculations of K'^^^. Finally, the 
excitation energy of the constrained isoscalar monopole vibration is computed as 



V -DM 

D. Comparison with the expressions in the nonrelativistic sum-rule approach 

We would like to point out that there is a parallelism between the scaling and constrained 
expressions for the GMR energy derived above for the relativistic model and the known 
results in the nonrelativistic formalism of sum rules. That is, one can define suitable average 
excitation energies of the resonance such as 



and E,^.^ (16) 

\ mi y m_i 

in terms of ratios of the integral moments = J^doj cu^SUj) of the RPA strength function 



S{uj). In the nonrelativistic RPA theory of the GMR [4J-|46|, it has been demonstrated that 
the cubic energy- weighted moment = J^duju^S{u) of the RPA strength function is 
exactly equivalent to the second derivative of the energy obtained from a monopole scaling 
transformation. It is also known that the RPA energy- weighted moment mi = J^duuS^u) 
is equivalent to 

mi = —A{r^) (17) 
m 

("Thouless theorem"), and that the exact RPA inverse energy-weighted moment m_i = 
du S (u) / u can be obtained from a constrained calculation of the Hartree-Fock ground 
state: 

("dielectric theorem"). Taking into account these results and the fact that m\ = {2/iv?) B'^j, 
one sees an interesting analogy between and Ei of the nonrelativistic RPA formalism of 
sum rules and -E^^^ of equations ([9]) and (E]) and E'^/^" of equations (ITSl) and (HM given 
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previously for the GMR average excitation energy in the relativistic model, though to our 
knowledge the sum-rule approach has not been derived in the relativistic RPA theory. 
It is worth mentioning that the i nteg ral moments of the nonrelativistic RPA strength 



function obey certain inequalities |44|-|46| . In particular, one has 



El < E, = — < ^3, (19) 
mo 

where Ex is the centroid energy of the GMR evaluated from the ratio of the mi moment to 
the mo moment. That is, Ei is a lower limit for the mean energy (centroid) of the resonance, 
whereas E^ is an upper limit. It is clear from the definition of the and m_i moments 
that E^ and Ei explore different energy domains of the collective excitation, i.e., E^ is more 
influenced by the high-energy region of the strength function S(u)), whereas Ei is more 
sensitive to the low-energy components of S{u). It can be shown [4^] that the quantity 



a = -^Ei-E! (20) 

provides an upper bound for the width of the giant monopole resonance in the RPA ap- 
proximation. The calculations carried out for stable nuclei show that E^ and Ei become 
increasingly close with increasing nuclear mass. Therefore, in a stable heavy nucleus the 
value of the jwidth a is small and the monopole strength becomes concentrated around a 
single peak [44j, which is the situation found in experiment. Following equation (120|) and 
the above comments, in the discussion of our results the spread between the calculated E]^/^^ 
and Ef}^ values will provide an indication about the width of the GMR. 



III. RESULTS AND DISCUSSIONS 

A. Comparison of relativistic ETF and Hartree-RPA results 

First, in this section we test with some examples the results of the relativistic extended 
Thomas-Fermi (RETF) approach vis-a-vis relativistic Hartree and RPA calculations. To 
this end, we report in table 1 semiclassical RETF results for the binding energy and the 
excitation energy of the GMR in the nuclei ^°*Pb and ^°Zr. We also include in the comparison 
the value of the neutron skin thickness ARnp = Rn — Rp (the difference between the neutron 
and proton rms radii of the nucleus) because this quantity is very sensitive to the density 
dependence of the symmetry e nerg y 0, 12, [t^, 71|. The results have been computed with 



the NL3 mean field interaction [72| supplemented by the nonlinear Ay coupling between the 
cj-meson and p-meson fields jsj. As pointed out in the Introduction, this model, through 
the isoscalar-isovector meson coupling, allows one to change the symmetry energy leaving 
unchanged the properties of the EOS of symmetric nuclear matter. By modifying the value 
of Ay one can investigate the changes in finite nuclei due to the variation of the density 
dependence of the symmetry energy almost without affecting the previous quality of the 
model for binding energies and charge radii jS-lO, 41, 43 



For the coupling constant Ay we use the values 0.00, 0.01, 0.02, and 0.03 (models denoted 
as NL3.00, NL3.01, NL3.02 and NL3.03 in table 1). It is known that the binding energy 
of nuclei constrains an average between the symmetry energy at saturation and the surface 
symmetry energy, rather than the individual value of the symmetry energy at saturation 
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TABLE I: Relativistic Hartree-RPA and relativistic extended Thomas-Fermi values obtained with 
the NL3 model plus the isoscalar-isovector coupling Ay for the binding energy per particle (MeV) , 
neutron skin thickness Ai?„p = i?„ — Rp (fm), and excitation energy of the GMR (MeV) in 208p|-, 
and ^'^Zr. The GMR centroid energies nii/niQ are taken from [41]. 



208pb. 

quantal Hartree-RPA semiclassical RETF 



Model 


B/A 




mi/niQ 


B/A 






TTiSCal 


NL3.00 


7.87 


0.28 


14.32 


7.99 


0.23 


13.92 


14.58 


NL3.01 


7.89 


0.25 


14.43 


8.01 


0.20 


13.97 


14.61 


NL3.02 


7.91 


0.22 


14.57 


8.02 


0.17 


14.07 


14.70 


NL3.03 


7.91 


0.20 


14.74 


8.03 


0.15 


14.20 


14.83 


90Zr: 




















quantal Hartree-RPA 




semiclassical RETF 




Model 


B/A 




mi/rriQ 


B/A 


ARnp 


ijicons 




NL3.00 


8.69 


0.11 


18.62 


8.91 


0.10 


19.04 


19.53 


NL3.01 


8.69 


0.10 


18.67 


8.91 


0.09 


19.05 


19.54 


NL3.02 


8.70 


0.09 


18.69 


8.92 


0.08 


19.08 


19.58 


NL3.03 


8.70 


0.08 


18.75 


8.92 


0.07 


19.12 


19.62 



(8l-[l2l|. Thus, as done in earlier literature [8-10, 41, 43 1, we constrain the symmetry energy 
S{p) of the models to take a value of 25.67 MeV at a Fermi momentum kp = 1-15 fm~^ 
(i.e., at a subsaturation density p ~ 0.10 fm~^). This implies that the coupling constant Qp 
of the NL3 interaction has to be suitably adjusted for each value of the isoscalar-isovector 
Ay coupling. As a consequence, the value S'(po) of the symmetry energy at the saturation 
density po changes with the Ay parameter. We find the results 5'(po) = 37.4, 35.0, 33.2 and 
31.7 MeV when we vary Ay from 0.00 to 0.03. That is, the symmetry energy S{p) of the 
model around saturation becomes softer (it increases more slowly with the nuclear density) 
when the value of the Ay parameter is larger. Further details on the behaviour of S{p) in 
the NL3 model supplemented by the Ay coupling can be found in (lo| . 

Table I shows that by increasing the Ay coupling the neutron skin thickness of a neutron- 
rich nucleus is significantly reduced, i.e., the models with a softer symmetry energy give rise 
to smaller neutron skins. Here, this is achieved while maintaining the proton rms radius 
of the nucleus essentially intact. For example, in the quantal Hartree calculation with the 
NL3.00 model we find Rp = 5.466 fm for ^'^®Pb and 4.199 fm for ^"Zr, while the respective 
values using NL3.03 are 5.478 fm and 4.210 fm. The changes in Rp are thus less than 0.5%. 
Moreover, one sees in Table I that also the binding energies are little altered by Ay. Both 
^"^Pb and ^°Zr are slightly more bound at large Ay but the changes in B/A are again no 
larger than 0.5%. Thus, the predictions for charge radii and binding energies of the original 
NL3 model are not seriously compromised by the additional Ay coupling, which otherwise 
largely softens the symmetry energy. 

From table 1 we observe that in the semiclassical RETF approach the binding energies of 
^°^Pb and ^"^Zr are larger and the neutron skins smaller, than in the corresponding quantal 
Hartree calculations. In principle, this fact is not to be considered a lack of accuracy of the 
Thomas-Fermi-like methods because semiclassical and quantal calculations must differ by 



10 



shell effects. It is well known that the semiclassical ground-state calculations by construction 
smooth out the quantal shell effects and provide the average part of the quantal binding 
energy |7i-[75|. Actually, in the semiclassical calculations one works with the particle densi- 



ties and completely avoids to calculate single-particle wave functions. By this same reason, 
the semiclassical approach yields smooth average densities and energies that do not show 
the quantum shell oscillations. The nuclear ETF functionals with gradient corrections up 
to order always produce some overbinding (between a ~2% in heavy nuclei and a ~10% 
in light nuclei) and yield smaller rms radii (between ~1% and ~4%) than the quantal mean 



field calculations [53|, |56|, |57|, |7^. The quantal corrections, which are more prominent in 
magic nuclei (our situation in table 1), can eventually be added perturbatively to the semi- 
classical binding energies using techniques such as the Strutinsky shell-correction method 



75l | or the expectation value method [56|, l57|, [74 



Our focus here is somehow more about the change of the results with the softness of 
the symmetry energy than on the absolute values themselves. In this respect, one has to 
note that the quantal and semiclassical results reported in table 1 show both practically 
the same variation as a function of the Ay coupling. The binding energies remain similarly 
constant and the neutron skins decrease by almost the same amount of 0.08 fm in ^osp]^ 
and of 0.03 fm in ^°Zr in both the Hartree and RETF calculations when Ay is changed 
from 0.00 to 0.03. Therefore, the RETF approach predicts practically the same variation 
with Ay as the quantal Hartree calculations. This implies that the rate of change of the 
results with the softness of the symmetry energy in the present model is well described by 
the RETF method. As a further test, which is motivated also because later we will study 
RETF results along isotopic chains, in figured] we plot for Pb isotopes the difference between 
the value of the neutron skin thickness of the nucleus calculated with Ay = 0.00 and with 
Ay = 0.03, as a function of the parameter / = (A^ — Z)/A. We observe by comparison 
with the corresponding quantal Hartree results, that the RETF approach performs quite 
satisfactorily along the Pb isotopic chain in predicting the right size of the modification 
induced by the change of the symmetry energy of the model. 

With regard to the average excitation energies of the GMR, table 1 shows the scaled and 
constrained values (i.e., -E^^' of section [HB] and E^^/^^ of Section HTCj) for ^osp]^ 902j~ 
calculated in the RETF approach. We compare our RETF results with the values of the 
centroid energy mi/mg of the GMR extracted from the relativistic RPA (RRPA) strength, 
which have been reported in [4l[. We see that for the ^osp]^ nucleus, the RRPA centroid 
energies very nicely lie in between the RETF constrained and scaled estimates, as it is 



required in the nonrelativistic sum-rule approach to the RPA |44j-|46l| (recall equation (fT9|) 
of section Hi Dl) . In the case of a lower- mass nucleus such as ^°Zr, the semiclassical excitation 
energies are not expected to be as successful as for the heavier ^"^^Pb. Indeed, we see in 
table 1 that the RETF constrained and scaled estimates in ^°Zr overestimate a little the 
RRPA centroid energies. The discrepancies are nevertheless moderate, with differences no 
larger than 2% for E'^^^^^ and 5% for E^^f^ with respect to the RRPA mi/rriQ value of ^"^Zr. 
And, as it happens in the case of ^°^Pb, the semiclassical average excitation energies of ^°Zr 
as a function of the Ay coupling display the tendency of the RRPA excitation energies. In 
summary, we can conclude from the results discussed in the present section that the RETF 
approach is a reasonable starting point to study the influence of the density dependence of 
the symmetry energy on the GMR of neutron-rich nuclei. 
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FIG. 1: Difference between the values of the neutron skin thickness ARnp obtained with Ay = 
and with Ay = 0.03 for lead isotopes. The results of the RETF method are compared with the 
results of quantal Hartree calculations. 

B. GMR in heavy neutron-rich nuclei 

As mentioned in the Introduction, we found in |47|] within the Skyrme-Hartree-Fock 
framework that the semiclassical calculations are able to describe realistically the tendencies 
of the excitation energy of the GMR in regions away from the stability valley, especially for 
heavy nuclei such as the lead isotopes where a semiclassical approach performs better. In 
the present section we analyze with the RETF approach to the relativistic mean field model 
the effects of the density dependence of the symmetry energy on the excitation energy of the 
GMR in heavy neutron-rich nuclei by contrasting the results for the nucleus ^^^Pb, which 
we take as an example of a heavy system with a large neutron excess away from stability, 
with the results for the stable isotope ^°^Pb. 

We display in figure |2]for ^osp]^ ^^^^ figure [3]for 266p]^ ^j-^g calculated values of the energy 
per particle E/A, the neutron skin thickness ARnp = Rn — Rp, the compression modulus of 
the finite nucleus Ka and the average excitation energy of the GMR as a function of the 
coupling constant Ay in the relativistic model. We recall that the coupling Ay induces a 
change of the density dependence of the symmetry energy. The model has a stiff symmetry 
energy at Ay = 0.00, while it has a soft symmetry energy at Ay = 0.03. The RMF 
Lagrangian supplemented by the nonlinear interaction between the u and p meson fields 
leaves the binding energy of the nucleus almost independent of Ay, at least for moderate 
values of this coupling. We note from the results for 266p]^ figure [3] that this happens 
not only in the case of stable nuclei, which we have seen previously for ^°^Pb and ^°Zr in 
table 1, but also in the case of much more neutron- rich nuclei. The neutron thickness ARnp 
decreases almost linearly with an increase of Ay in both ^''^Pb and ^^^Pb. The reduction of 
the ARnp value at Ay = 0.03 amounts to about a 30% of the original neutron skin thickness 
at Ay = 0.00. The proton rms radius Rp of ^osp^ ^^d 266p]^ j^^^g i^qq^ found to remain 
nearly constant with Ay, and therefore the decrease of ARnp is due almost exclusively to 
the decrease of Rn- It is worth mentioning that the alluded variation of the neutron skin 
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FIG. 2: Dependence with the value of the isoscalar-isovector coupUng Ay of the energy per particle 
(a), neutron skin thickness (b), finite nucleus compression modulus (c) and average excitation 
energy of the GMR (d) for the ^'^^Pb nucleus. 



thickness of ^°^Pb in the present calculations is close to roughly a 1% change in the neutron 
rms radius i?„ of this nucleus. That is, the range of variation we have allowed in the density 
dependence of the symmetry energy by changing Ay between 0.00 and 0.03 induces a change 
in Rn of 208pi3 around 1%. In this sense, it is consistent with the uncertainty that will remain 
in this observable after the PREX experiment at the Jefferson Laboratory [76|], which aims 
to measure Rn of ^osp]^ 2% accuracy via parity-violating electron scattering. 

The average excitation energy of the GMR calculated in the RETF approach involves the 
knowledge of the finite nucleus compression modulus K^, or restoring force of the resonance. 
We display in figures [2]^c) andlU^c) for ^osp]^ g^j^^ ^^^Pb, respectively. Let us discuss first 
the change of Ka with the mass number A in isotopes. We see from the results for ^°^Pb 
and 266p]^ that the finite nucleus compression modulus decreases when A increases. This 
behaviour may be easily understood from the so-called leptodermous expansion of Ka (77j : 



Ka = Ko + KsurfA"'/' + Kri' + KcouiZ^A-^/' + ■ ■ ■ , (21) 

where Kq is the compression modulus of symmetric infinite nuclear matter and -ft'surf , -^r and 
Kcoui are the leading surface, volume-symmetry and Coulomb corrections. The coefficients 
i^sii rf , Kr, and -ft'coui in the scaling model are negative, and in particular K^. may be large 



33l . l63l . |77H80| . Within an isotopic chain, the nuclear charge Z is fixed and therefore the 
neutron-proton asymmetry / = (A^ — Z)/A increases when A increases. Then, although the 
surface and Coulomb contributions in fl2Tl) become a little less negative with increasing A, 
the increase of the negative volume-symmetry contribution Krl^ dominates. Therefore, Ka 
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FIG. 3: Dependence with the value of the isoscalar-isovector coupUng Ay of the energy per particle 
(a), neutron skin thickness (b), finite nucleus compression modulus (c) and average excitation 
energy of the GMR (d) for the ^^^Pb nucleus. 



in isotopes decreases with larger A. We stress that the expansion fl2Tl) is provided here for 
indicative purposes only and that in all results we have computed Ka self-consistently. 

The scaling value K^^^^ of the finite nucleus incompressibility, given by equation (|8]) of 
section IIIBt is displayed in figures |2t^c) and |3t^c) by the upper curves with solid symbols. For 
both isotopes ^osp]^ g^j^^ ^^^Pb, K^^^ shows an upward trend with larger Ay- This variation 
of K^^^ with Ay is due to the softening of the density dependence of the symmetry energy. 
It modifies the value of K^. and consequently Ka (of course, K^- does not describe alone 
the whole effect as there are surface-symmetry and other corrections neglected in fl?I]) [63|). 
The present models have all Kq = 271.5 MeV, while the Kr coefficient changes from —699 
MeV in NL3.00 to —381 MeV in NL3.03. The recent empirical extractions of Kr suggest Kr 
values in this range. For example, isospin diffusion in nuclear reactions favours the constraint 
Kr = —500 ± 50 MeV the breathing mode of Sn isotopes suggests Kr = —550 ± 100 
MeV [2^] (or -500 ±50 MeV [79|) and neutron skins of nuclei point to Kr = -500l}^^ MeV 

The excitation energy -E'^'}^' of the monopole oscillation evaluated in the scaling approach 
is determined by the ratio of the restoring force K^^^ to the mass parameter Bm (section [ITll. 
As indicated, to replace Bm by its nonrelativistic limit B^ is a very good approximation. 
In practice, this implies that Bm varies with the parameters of the nuclear interaction 
essentially in the same way as (r^) of the nucleus does, see equation ffTTl) . When Ay is 
increased, (r^) of the matter distribution of the nucleus decreases because of the reduction 
of the neutron rms radius originated by the softer symmetry energy. As a consequence. 
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the mass parameter of the GMR (analogously, the mi sum rule (fT71) in the nonrelativistic 
language) also decreases. This effect combined with the growing tendency of K^'^^ with 
Ay implies that the scaling excitation energy of the GMR will be enhanced with a softer 
symmetry energy, which is confirmed in figures[2]^d) and[3](d) where E]^'^^ is plotted. We recall 
that this enhancement of E'^f'^ is due to the different density dependence of the symmetry 
energy alone because the incompressibility Kq of symmetric matter is forced to be the same 
in all the interactions considered here. 

Let us now analyze the results obtained through the constrained calculations described 
in section III C[ For the stable nucleus ^"^^Pb, we see in figure [2] that K'^^^ and Ef}^^ have 
a similar behaviour to the discussed scaling results as a function of the density dependence 
of the symmetry energy. The separation between E^j^^^ and E'^^^ is indicative of the width 
of the GMR resonance (see equation fl20|) ). and therefore this quantity is basically constant 
with the symmetry energy in ^°^Pb. In contrast to this situation in the stable ^osp]^ nucleus, 
in the more neutron-rich isotope ^^®Pb we see in figure [3] that the change of K'^^^ and 
El""^^ with the symmetry energy exhibits a different pattern from the corresponding scaling 
results. For ^®^Pb, both quantities K'^^^ and -E^/"*^ '^'^^ only do not increase with Ay, but 
they decrease significantly. 

To interpret the reduction of E'}^^^ and K'^^^ with Ay in ^^^Pb, it will be helpful to 
rewrite the expression Eff^^ = {AK'^^^ / BmY^"^ of equation ffTSj) by approximating Bm with 
its nonrelativistic form = mA (r^), and to write K'^^^^ using the first expression given 
in equation (IMI) . This results in 



Em-^J ^J^, (22) 



j^(3ii?,/Sr,),.o 



which in terms of the defined quantities mi and m_i has the same form that E^^ takes 
in the nonrelativistic RPA formalism of sum rules (i.e., Ei of equation ( I16p ) |44j-|46| . It is 
to be noted that the denominator of is very sensitive to the presence of loosely bound 



nucleons in the nucleus 47|], which is precisely the situation that occurs in ^^^Pb. When 
one is approaching the neutron drip line, the negative chemical potential /i„ of neutrons 
becomes close to vanishing, implying that the system contains weakly bound neutrons. 
These neutrons are very soft against compression and in the constrained calculations they 
give rise to a large increase of the absolute value of {dR'^/dri)rj=Q, i.e., of m_i in fl22l) . This 



effect has been found in [47|] in constrained HF (and ETF) calculations for nonrelativistic 
Skyrme forces, where it can be related to the enhancement of the RPA strength in the low- 
energy region produced by the increase of the transitions to the continuum at large neutron 
numbers. The same effect is found in the recently published QRPA and CHFB Skyrme 
calculations, which take account of the pairing correlations, performed in the light isotopic 
chains of Ca and Ni [37|. In our present calculations, the increase of m_i at large neutron 
numbers is accentuated the softer is the symmetry energy (note that parameter sets ha ving 
a softer symmetry energy have also weaker binding energy of the last bound neutrons 10 
and therefore reinforce the effect). This explains the decreasing trend of E'^^ with Ay in 
266p]^ noticed in figure |3l The decrease of K'^'^'^ with Ay seen in the same figure is similarly 
understood by recalling equation ( IT^ that relates K'^'^^ with {dR^ dri)~}^Q. 
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C. GMR along the Pb and Zr isotopic chains 



In this section we would Uke to investigate at which values of the neutron excess of the 
nucleus the effects coming from the symmetry energy become more relevant in the excitation 
energies of the GMR in isotopic chains of heavy (Pb) and medium (Zr) mass. Within the 
Skyrme-Hartree-Fock framework, in 47| we analyzed the excitation energy and the resonance 
width of the GMR in some isotopic chains with the sum-rule approach (i.e., through scaling 
and constrained calculations). In the case of a small atomic number such as the Ca isotopes, 
the semiclassical ETF results were able to describe the global variation of the quantal HP 
excitation energies but could not reproduce them quantitatively. However, the ETP results 
for the GMR excitation energies become progressively more accurate with larger atomic 
number. In the case of the Pb chain, the ETP scaling and constrained excitation energies 
were found to be in very good agreement with the quantal HP values on the quantitative level 
from the proton to the neutron drip line (the differences between ETP and quantal values 
not exceeding 0.4 MeV) Thus, from the experience in the nonrelativistic calculations, as 
well as suggested by the comparison of RETP and RRPA results made in table 1, the present 
semiclassical predictions for the GMR are expected to be sufficiently realistic, and rather 
accurate quantitatively (compared to quantal calculations) in the case of the Pb isotopes. 

We depict in figure H] the scaling and constrained average excitation energies of the GMR 
for the Pb isotopic chain. Here we only use the values 0.00 and 0.03 for the nonlinear u-p 
interaction Ay, representing stiff and soft density dependences of the symmetry energy. In 
the semiclassical calculations, the vanishing of the neutron chemical potential pn defines the 
location of the neutron drip line. We see in figure S] that the latter is reached at a smaller 
mass number when the model has a soft symmetry energy, in agreement with the earlier 
observations of Note that shell effects are absent in these semiclassical predictions of 
the drip line (a quantal Hartree calculation with NL3 predicts the Pb neutron drip line 
a.t A = 266). The proton drip line occurs at neutron number not much larger than Z. 
Therefore, because one deals with almost symmetric matter, the position of the proton drip 
line is unaffected by the density dependence of the symmetry energy. 

The calculated GMR excitation energies show in figure H] a downward tendency with 
increasing mass number A when one proceeds along the isotopic chain, by similar reasons 
to those pointed out in our previous discussion about the ^osp]^ ^^-^^ ^^^Ph isotopes. The 
downward trend is more noticeable for E'^^^ than for E'^'}'^^ in the isotopes with larger neutron 
numbers. The scaling excitation energies E]^'^^ obtained using a parameter set with a soft 
symmetry energy (Ay = 0.03) are systematically higher in the whole isotopic chain than if 
the symmetry energy is stiff (Ay = 0.00). Similarly to the scaling results discussed for ^ospj-^ 
and ^^^Pb, this behaviour of Ef^^^ is a consequence of the fact that the matter radius of the 
nucleus is smaller and the scaling incompressibility K^^^ larger for a soft symmetry energy. 

At variance with the uniformity displayed by the E]^'^^ results, the constrained estimate 
E'^'^ of the GMR excitation energy has a distinct feature with the symmetry energy below 
and above A ~ 254. If the mass number is below A ^ 254, E^j"^^ is larger when the parameter 
set has a soft symmetry energy. However, this tendency is reversed from A ~ 254 to the 
neutron drip line when the isotopes are more neutron rich. Actually, in a given nucleus 
the Ay parameter reduces its neutron rms radius, which means that the system becomes 
on average more compact. As long as the nucleus does not be very neutron rich, if it is 
more compact, the effect of the constraining parameter t] on the nuclear rms radius is less, 
and therefore the absolute value of {dR^/dvi)j^=Q in the denominator of equation fl22l) for 
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FIG. 4: Mass-number dependence of the excitation energy of the giant monopole resonance in the 
Pb isotopic chain for the NL3 model with Ay = (stiff symmetry energy) and Ay = 0.03 (soft 
symmetry energy), in the scaling and constrained approaches. 



E'^'^ is smaller. This explains the larger magnitude of E'^"'^^ at Ay = 0.03, compared to 
Ay = 0.00, in the Pb isotopes lighter than A ~ 254. But a point is reached in mass number 
(A ~ 254 here) where the isotopes are so neutron rich that the last neutrons become very 
weakly bound. These outer neutrons are very soft against compression. As a result of their 
specific contribution to when one applies an external constraint, the absolute value of 
{dR'^/dr])fj=o increases sharply instead of decreasing. This originates the large reduction of 
Elf^^ seen in figure H] for the more neutron-rich isotopes. The reduction of E'^^ at large 
neutron numbers is enhanced when the symmetry energy is soft compared to the case of a 
stiff symmetry energy because the outer neutrons are more weakly bound in the former case. 
A comparable decrease of the excitation energy of the GMR, however, is not present in the 
scaling calculations because the discussed effect does not arise in a self-similar scaling of the 
density of the nucleus. As we mentioned in section III D\ the scaling and constrained average 
energies pertain to different energy regimes of the GMR. While E]^'^^ informs more about 
the relatively high-energy contributions to the RPA strength, -E^^"*^ carries more information 
about the low-energy sector which is the region that shows the more noticeable changes in 
largely neutron-rich nuclei because of the presence of weakly bound neutrons. 

We may take a = |[(^17')^ - (^m''')T^^ as an estimate of the width of the GMR, in 
analogy with the nonrelativistic framework 44-46[ (see section Hi DP . From the separation 
observed in figure H] between E^j^'^^ and E^^^^, at either Ay = 0.00 or Ay = 0.03, the width of 
the resonance is expected to be almost independent of the mass number and of the density 
dependence of the symmetry energy for all of the Pb nuclei located in the region between 
the proton-rich side of the isotopic chain till A ~ 230-240. However, E^^^^ starts to show a 
pronounced departure from E^^f"^ after A ~ 240 and, as a consequence, the resonance width 
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FIG. 5: Mass-number dependence of the excitation energy of the giant monopole resonance in the 
Zr isotopic chain for the NL3 model with Ay = (stiff symmetry energy) and Ky = 0.03 (soft 
symmetry energy), in the scaling and constrained approaches. 



should display a sharp increase for the Pb isotopes with large neutron numbers. We see 
from figure H] that this effect is predicted to be stronger when the nuclear interaction has 
a soft symmetry energy. As discussed in [47] in nonrelativistic calculations with Skyrme 
forces, the value of the estimate a of the GMR width provides a qualitative idea about the 
distribution of the RPA strength. A small width a indicates that the RPA strength of the 
GMR is basically concentrated in a narrow single peak, while a large width a suggests that 
the peak is broad, or even that the RPA strength may be fragmented into several peaks. The 
large value of m_i oc —{dR'^/dr])rf=o (see equation (|22l) ) that one finds in the constrained 
calculations for the more neutron-rich Pb isotopes, points towards an enhancement of the 
RPA strength in the low-energy region in these nuclei, due to transitions from the last 
weakly bound single-particle levels to the continuum [43|. (This effect, which has a quantal 
origin, is accounted at least on average by the semiclassical calculations [471].) The discussed 
scenario is actually confirmed by the QRPA calculations of [s^] where the fragmentation as 
well as the enhancement of the low-energy part of the RPA strength in Ca and Ni isotopes 
approaching the neutron drip line can be seen in figures 3 and 6. 

Lastly, we present in figure |5]the calculated results for the excitation energies of the GMR 
in the Zr isotopic chain. The variation of the excitation energies obtained in the scaling 
scheme as a function of the mass number A and of the Ay coupling is qualitatively similar 
to the pattern exhibited by the Pb isotopes. That is, the excitation energies -E^^^ show a 
monotonic decrease with increasing mass number, and they are larger when calculated with a 
soft symmetry energy. The relative differences between the values of E^^f^ at Ay = 0.03 and 
at Ay = 0.00 are smaller compared to the chain of the Pb isotopes because of the smaller 
isospin content of the Zr isotopes. Concerning the results for the constrained excitation 
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energies in Zr, the value of E'^^ shows a stronger dependence on the mass number A at 
large neutron numbers than the scaling value -E^*^'. As long as A < 120, the excitation 
energy E'^^^ for zirconium is to a large extent independent of the symmetry energy. But 
afterwards {A > 120), E'^'^'^ becomes visibly lower for the soft symmetry energy. 

In consonance with the behaviour of Elf'^ and E^^^^ found in figure it is predicted that 
the width a of the GMR for the Zr nuclei should remain practically constant till A ~ 100. 
The width a would start increasing moderately between A ~ 100 and A ~ 120, because of 
the slight gradual separation between and E'^^^ in these isotopes. Finally, the resonance 
width would be much enhanced in the region of the most neutron-rich Zr isotopes, where 
the effect is largely reinforced by having a soft symmetry energy in the nuclear interaction. 



IV. SUMMARY AND CONCLUSIONS 

Though the excitation energy of the isoscalar giant monopole resonance is mainly ruled 
by the compression modulus of symmetric nuclear matter, it is also affected by the symmetry 
energy in neutron-rich nuclei. In this paper we have analyzed the influence of the density 
dependence of the symmetry energy on the average excitation energy of the GMR. In order 
to do this study, we have used a relativistic mean field model supplemented by an additional 
u-p nonlinear interaction driven by a coupling constant Ay- This additional coupling allows 
one to generate different parameter sets that preserve the properties of the symmetric matter 
and differ among them, only, in the density dependence of the sjTiimetry energy, which is 
softened by increasing Ay. As a function of the coupling constant Ay, the model predicts 
in finite nuclei almost the same binding energy and proton rms radius as in the case where 
Ay = 0, while the neutron rms radius gets significantly reduced. 

The excitation energy of the isoscalar giant monopole resonance has been computed 
using the semiclassical relativistic extended Thomas-Fermi theory. We obtain the energy 
of the breathing mode in two different ways: through the scaling model and performing 
constrained calculations. It is known that these calculations with nonrelativistic Skyrme 
forces nicely reproduce the average energies E^ and Ei obtained from the more microscopic 
RPA calculations, not only for stable nuclei but also for largely neutron-rich nuclei |47|. 
In our study, we have calculated the scaling and constrained semiclassical average energies 
with the relativistic model along the isotopic chains of Pb and Zr. We can also obtain some 
qualitative information about the distribution of the RPA strength from the resonance width 
computed using these two semiclassical estimates of the excitation energy. 

As it is known, the excitation energy of the breathing mode decreases, at least for stable 
nuclei, following roughly an A~^^^ law |77|] due to the geometrical enhancement of the mass 
denominator. However, when the neutron-proton asymmetry increases, departures from this 
behaviour appear. The reason is that the isospin-dependent part of the finite nucleus in- 
compressibility Ka plays a non-negligible role by reducing the restoring force and making 
nuclei with large neutron numbers softer than the stable isotopes. This effect in the aver- 
age excitation energy is present in the scaling estimate but it is particularly strong in the 
constrained calculation. When the neutron-proton asymmetry increases, the outer neutrons 
become more and more weakly bound starting to behave almost as extremely asymmetric 
nuclear matter, which is very soft compared to symmetric matter 8ll|. This effect makes 



the compression modulus K'^'^'^ smaller and, therefore, the average excitation energy of the 
GMR decreases. As discussed in the case of constrained HF calculations with Skyrme forces 



47( 1 . this scenario reflects the important contribution to the low-energy part of the RPA 
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strength distribution coming from the particle-hole transitions to the continuum from the 
weakly bound neutron energy levels in nuclei having large neutron numbers. A similar sit- 
uation has been found recently by Capelli, Colo and Li 37| in QRPA and constrained HFB 
Skyrme calculations performed along the light isotopic chains of Ca and Ni. 

For a given nucleus the excitation energy of the GMR also depends on the density content 
of the symmetry energy of the model used to compute it. In stable nuclei the scaled and 
constrained estimates of the excitation energy are larger for parameter sets that have a 
symmetry energy with a soft density dependence. In the case of the constrained estimate, 
this tendency is reversed in the nuclei with large neutron-proton asymmetry. The reason 
lies in the fact that in these nuclei the outermost and weakly bound neutrons are very 
compressible and consequently the constrained compression modulus K'^'^^ becomes lower 
when the density dependence of the symmetry energy is soft (the outer neutrons are then less 
bound) than when the symmetry energy is stiff. Thus, the average excitation energies based 
on the inverse energy-weighted sum rule are predicted to be significantly reduced for large 
neutron numbers when computed with a parameter set that has a softer density dependence 
of the symmetry energy. 

In conclusion, the density dependence of the symmetry energy has a moderate influence 
on the value of the excitation energy of the isoscalar giant monopole resonance in stable 
nuclei and in nuclei that are not far from stability. It has a larger impact on the GMR 
of isotopes with large isospins, where the calculations reveal that the constrained energies 
become much reduced in models having a soft symmetry energy. The results point towards a 
considerable enhancement of the low-energy part of the RPA strength and to likely different 
structures in the excitation spectra of the GMR in unstable neutron-rich nuclei if the nuclear 
symmetry energy is soft. It may be worth pursuing a detailed characterization of these 
predictions related with the softness of the symmetry energy by computing the strength 
distribution in a few selected isotopes at large neutron numbers with the more demanding 
but more microscopic QRPA theory of the GMR response. Unfortunately, it is difficult 
that the indicated effects be tested in the laboratory as they become prominent at neutron 
numbers beyond present reach. Such nuclear systems are otherwise expected to exist in 
astrophysical conditions, such as those prevailing in the inner crust of neutron stars, where 
exotic nuclear clusters can be stabilized by an enveloping gas of neutrons and electrons (50| . 
It is thought that low-energy collective excitations play an important role in astrophysical 
processes that involve nuclei far from stability [82[, and in that context the modification of 
the low-energy region of the strength distribution of the monopole resonance by the softness 
of the symmetry energy could have implications. 
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Appendix A: Equations of the relativistic extended Thomas-Fermi model 



We have introduced the energy density of the relativistic nuclear mean field model supple- 
mented by an isoscalar-isovector mixed u-p interaction in equation ([1]) of the text. Following 
, in the relativistic extended Thomas-Fermi approach the expression of the nucleonic 



part £ is written as £ = £q + £2-, where 



y- 



g=n,p 



4 , k-pq + epg 



rrr 



(Al) 



and 



^2=Y. [C'i,(fcF„m*)(Vp,)2 + C2,(A;F„m*)(Vp,-Vm*) + C3,(A;F„m*)(Vm*)2]. (A2) 



i}=n,p 



The leading term £q is the contribution of the pure Thomas- Fermi part, while the term £2 
contains the gradient corrections of order fi^ . The variables /cpg, and ep, are, respectively, 
the local Fermi momentum fcpg = (Svr^Pg)^/^, the Dirac effective mass m* = m — $, and the 
dispersion relation epg = -|- m*^)^/^ of the effective particle, with g = n for neutrons and 
g = p for protons. The coefficients Cjg stand for the following analytic functions of /cpg and 
m* = m — $: 



C 
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2q 



94^3 ,2 



m 



Qk-pqe^ 



In 



Fg 



epq + 2fcFg In 



k^q + epg 



k^q + CPg 



3g 



'^Pg 



CPg 
k^q 



Fq \ fepg + epg 



'*'Pg 



(A3) 



The RETF ground-state densities and meson fields are calculated by solving by numer- 
ical iteration the Euler-Lagrange equations associated to the energy density ([T]), with the 
constraint of baryon number conservation (which is imposed through Lagrange multipliers 
Pq, the chemical potentials of neutrons and protons): 

1 



dC2q dCsg 



dm* 



dpq 



Vm*)^ = m + pq. 



(A4) 
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(A-m2)$ = -gUp 



, ^$2 _ A$3 
' 2! 3! 



(A - ml)W = ~gl (p - 2kyB^W) , 
{A-ml)B = -gl{p,-2AvBW^) 
AA = -Pp. 



(A5) 

(A6) 
(A7) 
(A8) 



We have r^g = +1 for protons and = —1 for neutrons. The baryon density is p = Pp + pn, 
while ps = |(pp — Pn) is the isovector density of the nucleus. Finally, the RETF expression 
of the scalar density ps entering equation f lA5l) is given by 



ps 



5£ S£o ^ 5£2 ^ 

6m* 6m* 6m* ^'^^ 



Em* 
2^ 



1 *2 1 k-pq + e-pq 

KFq^Fq — m lu 



-E 



C2qApq + 2C3,Am* + 



dC2q dC 



dpq dm* 



(Vp,) 



da 



+ 2^(Vp,-Vm*) + 

dpq 



da 



39, 



(9m* 



Vm 



*n2 



(A9) 



Appendix B: Klein-Gordon equations of the scaled meson fields 



In this appendix we illustrate how to calculate the derivatives of the scaled meson fields 
upon the scaling parameter a that we have used to represent the collective coordinate of 
the monopole oscillation of the nucleus in the scaling approach. In view of equation (E]) for 
the scaled baryon density p^, i.e., Pa{,f) = a^p{ar), and of the field equation (1A6P obeyed 
by the omega-meson field W, the scaled field Wa fulfills the Klein-Gordon equation 



OA 

gUapiu)-^B^iuyW^iu] 



(Bl) 



where we have introduced the coordinate u 
respect to the scaling parameter a we have 



= ar. On differentiating equation fIBip with 



A„ - 



da 



-9v[P + 



2ml 4A 



gl 



+ 



a 



4A, 



■B^W, 



dB^ 
da 



2Av ^2 9Wa 

a"^ " da 



(B2) 



If we now set a = 1, the solution of this equation allows us to obtain the value of the 
derivative dWa/ da\a=i- One can work out the equations for the scaled isovector p- meson 
field in the same way. In practice, one simply has to replace {Wa, Ba, p, m„, g^) by 
{Ba, Wa, Ps, rup, gp) in equations (IBip and (IB2p . 
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In the case of the scalar field one has to take account of the terms that arise due to the 
fact that the scalar density ps is itself a function of the scalar field. Following the same steps 
as above, from the Klein-Gordon equation 

(a„ - ^) $„(n) = -aglpfiu) (B3) 

for the scaled scalar field one readily arrives at the equation 

whose solution at a = 1 provides the value of d^a/da\a=i. 
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